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There is a general belief, reinforced by statements in standard textbooks, that: (i) one can obtain 
the full non-linear Einstein's theory of gravity by coupling a massless, spin-2 field h a b self-consistently 
to the total energy momentum tensor, including its own; (ii) this procedure is unique and leads to 
Einstein-Hilbert action and (iii) it only uses standard concepts in Lorentz invariant field theory and 
does not involve any geometrical assumptions. After providing several reasons why such beliefs are 
suspect — and critically re-examining several previous attempts — we provide a detailed analysis 
aimed at clarifying the situation. First, we prove that it is impossible to obtain the Einstein-Hilbert 
(EH) action, starting from the standard action for gravitons in linear theory and iterating repeatedly. 
This result follows from the fact that EH action has a part (viz. the surface term arising from second 
derivatives of the metric tensor) which is non-analytic in the coupling constant, when expanded in 
terms of the graviton field. Thus, at best, one can only hope to obtain the remaining, quadratic, 
part of the EH Lagrangian (viz. the F 2 lagrangian) if no additional assumptions are made. Second, 
we use the Taylor series expansion of the action for Einstein's theory, to identify the tensor S ab , to 
which the graviton field h a b couples to the lowest order (through a term of the form S ab h a b in the 
(jj lagrangian). We show that the second rank tensor S ab is not the conventional energy momentum 

jy} ■ tensor T ab of the graviton and provide an explanation for this feature. Third, we construct the full 

nonlinear Einstein's theory with the source being spin-0 field, spin-1 field or relativistic particles 
by explicitly coupling the spin-2 field to this second rank tensor S ab order by order and summing 
up the infinite series. Finally, we construct the theory obtained by self consistently coupling h a b to 
the conventional energy momentum tensor T ab order by order and show that this does not lead to 
Einstein's theory. The implications are discussed. 
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A. Conventional wisdom 
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The two classical fields — electromagnetism and gravity — are described by a vector field and second rank symmetric 
tensor field, respectively. Considerations based on Lorentz group suggest interpreting them (when suitable restrictions 
are imposed) as corresponding to massless spin-1 and spin-2 fields. The vector field A4 couples to a conserved current 

^Jy Ji but does not contribute to this current. (That is, the photon does not carry charge.) In contrast, the tensor field 
is believed to be coupled to the energy momentum tensor; since the field itself carries energy, it has to couple to itself 
in a non linear fashion. (The situation is similar to Yang-Mills fields which carry isotopic charge and hence are non 
k> , linear.) It may, therefore, be possible to obtain a correct theory for gravity by starting with a massless spin-2 field 

5.^ 1 h a i, coupled to the energy momentum tensor T a j, of other matter sources to the lowest order, introducing self-coupling 
of h a b to its own energy momentum tensor at the next order and iterating the process. This will lead to a completely 
field theoretic description of gravity in a Minkowski background and is conceptually quite attractive. 

This attempt has a long history. The field equation for a free massless spin-2 field was originally obtained by 
Fierz and Pauli lj. The first attempt to study the consequences of coupling this field to its own energy momentum 
tensor seems to have been by Kraichnan in unpublished work done in 1946-47. The first published attempt to derive 
the non linear coupling is by Gupta and Kraichnan published some of his results soon after Feynman has 
provided a derivation [4| in his Caltech lectures on gravitation during 1962-63. The problem was re-addressed by a 
clever technique by Deser JjJ. (This problem and related ideas have been explored from several other points of view 
in literature; see e.g., [jj. We shall not discuss these approaches.) Virtually all these approaches claim to obtain not 
only Einstein's field equations but also the Einstein-Hilbert action. 



B And why it is suspect 

This result is widely quoted in literature (see, e.g., page 424 of 0) and, at first sight, seems eminently reasonable. 
However, deeper examination raises several disturbing questions, if the above result is really valid. 

• In the conventional derivations, the final metric arises as g a b = T]ab + A/i a 6 where A oc \ G has the dimension of 
length and h a f, has the correct dimension of (length) -1 in natural units with K = c = 1. The iteration is in powers 
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of A, starting from the zeroth order lagrangian Lq ~ (dh) 2 for spin-2 field, which has the dimension of (length) -4 . 
(We have dropped the tensor indices to simplify the notation.) The final result in all the published works is the 
Einstein-Hilbert lagrangian Leh = (l/4A 2 )i?. Since the scalar curvature has the structure R ~ (dg) 2 + d 2 g, 
substitution of g a b = rj a b + Xh a b gives to the lowest order: 

L EH oc -^R ~ (dh) 2 + \d 2 h (1) 
X A X 

Thus the full Einstein-Hilbert lagrangian is non-analytic in A! It will be quite surprising if, starting from (dh) 2 
and doing a honest iteration on A, one can obtain a piece which is non-analytic in A. At best, one can hope to 
get the quadratic part of Leh which gives rise to the T 2 action but not the four-divergence term involving d g. 

• To carry out this programme, one need to identify the energy momentum tensor T^ b for the graviton field h a b 
order by order in the coupling constant. At this stage, we are working in flat spacetime, Cartesian coordinates 
[with metric r\ a b = dia(— 1, 1, 1, 1)] with Lorentz group as the invariance group. If we are honest (and do not 
use anything we learnt in our general relativity course!), we must provide a prescription to find T^ b within this 
context. There is indeed a natural conserved second rank tensor which arises from Lorentz symmetry, usually 
called the canonical energy momentum tensor. This tensor, unfortunately, is not symmetric. It can be made 
symmetric but the procedure is not unique. For every choice of Tf? one can one obtain a nonlinear theory clearly 
showing that further choices are to be made somewhere along the line. 

• A sharper way of stating the above difficulty is the following: The same textbooks which assert that Einstein's 
theory can be obtained by coupling h a b to itself self consistently will also state in some other section (in ref. 
Q, this happens in page 467) that gravitational field does not have a well defined energy momentum tensor!. 
It will be rather strange if a unique energy momentum tensor exists for gravitational field order by order in 
the perturbation series but somehow "disappears" when all the terms are summed up. (The non-uniqueness of 
energy momentum tensor for Einstein's theory is well known and is extensively discussed in literature; see e.g., 
8J.) 

• In implementing this program, one needs to be clear whether general covariance is an assumption or a result. 
The starting point — Lorentz invariant field theory in flat spacetime with metric r/ a b — has no notion of general 
covariance. If the source of the final equations is an energy momentum tensor which is assumed to be generally 
covariant, it is equivalent to assuming that the left hand side of the equations is generally covariant. It is then no 
big deal to obtain Einstein's theory, if we are prepared to assume general covariance. [It is sometimes claimed 
that the gauge invariance of spin-2 field under h a b(x) — > h a b(x) + d a £t,(x) + db£, a (x), "becomes" the general 
covariance of the full theory. This is simply wrong; see the discussion around Eq. (|32H in Section III below.] 

• A term in the lagrangian proportional to Xh a bT ab where T ab is due to external matter fields (assumed to be 
independent of h a b to this order), will lead to the equation of motion of the type d 2 h = XT. Hence a coupling of 
the type Xh a bT ab is equivalent to requiring the source being T a b- Consider now the coupling of gravity to itself 
through a term of the type Xh a bS ab (h) where S ab explicitly depends on the graviton field h a b. When this term 
is varied with respect to h a b to get the equations of motion we will obtain two terms: S ab + (dS 1 ^ /dh a b)hij both 
of which will act as a source to gravity at next order. If we want the source to be the energy momentum tensor 
of graviton field, T a b, then the coupling cannot be of the form h a bT ab (h) since this will lead to the wrong source. 
Thus we need to find out the form of the tensor S ab — a question which does not seem to have attracted any 
attention in the literature. (We will see that S ab is an interesting object in its own right.) 

None of the previous derivations addresses these issues and most of them downplays the role of assuming general 
covariance. All these attempts make different tacit assumptions and it is difficult to judge which of these derivations 
can be thought of as "from first principles" in the sense that it is completely independent of our knowledge of the end 
result. This difficulty becomes apparent when one follows the details of many of these derivations. The technology 
used is very strongly influenced by the known final result. For example, Kraichnan's pioneering work explicitly uses 
a term like r) a bR ab (r]) (where rj a b is the Minkowski metric and R ab is the Ricci tensor) cleverly to obtain the result, 
[see equations 13-17 of |3j] in spite of the fact that R ab (r]) vanishes for the flat metric r/ a b ! It is impossible (at least 
for the author) to imagine that someone could have "guessed" this form for the action without knowing the result. 
Feynman's derivation also suffers from several shortcomings. To begin with, it is considerably less general than the 
one by Kraichnan since Feynman assumes a particular form for the matter action and a coupling. But more relevant 
to our discussion is the manner in which he constructs the solutions to a consistency condition (see the discussion in 
section 6.3 and 6.4 of Q). Since this approach assumes general covariance ( in the form X7 a T ab = 0) and relies heavily 
on constructing generally covariant scalars, it is predestined to give Einstein-Hilbert action. The by far cleverest 
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mathematical procedure was the one employed by Deser, in which he exploits the fact that, with a suitable choice 
of variables, the gravitational action becomes a cubic polynomial allowing the iteration to stop at a finite order. To 
achieve this mathematical economy, he has to start with the Palatini variational form (see his equation 2) based on 
the Lagrangian f ab R a b where f ab = \ / —gg ab is the preferred variable. It is no surprise that he obtains ^/—gR as 
the final result. (This is the only previous work that actually attempts an iteration; we shall comment on Deser's 
derivation in more detail in Section VI C.) 

In particular, the fact that the action for the final theory contains the second derivatives of the field is always put 
in by hand. Kraichnan's work has this explicitly; Feynman's derivation assumes a condition equivalent to general 
covariance to obtain solutions to a functional constraint and he explicitly chooses the "simplest" generally covariant 
scalar thereby getting y/—gR; Deser starts with an expression which is numerically same as y/—gR but since he uses 
the first order form of the variation, the question of second derivatives is not directly applicable until, of course, when 
the final result is obtained. Thus, we really do not know whether the original program (of coupling the field to its 
own energy momentum tensor order by order and iterating the result) can be made to yield any other form of action 
principle even when the field equations are the same. 

One may be tempted to argue that most of the issues and objections raised above are irrelevant in the strictly 
classical context. In classical general relativity, one could argue, what matters is the equations of motion and not 
action functional. This, however, is a rather restricted point of view and one needs to realise that the true world is 
quantum mechanical and if one can gain insight into the nature of quantum theory from the structure of classical 
action functional, it is worth exploring. Of course, quantum theory has taught us that action functionals are as 
important (if not more) as the field equations. In the case of gravity, there are two action functionals which are of 
primary relevance. The first is the Einstein-Hilbert action which uses the Lagrangian R^/—g and the second is the T 2 
action involving only the squares of the first derivatives. It can be shown that these two actions can be thought of as 
providing the momentum representation and the coordinate representation (respectively) of the theory and differ by 
a surface term which is directly related to the entropy of horizons in the semi classical theory and has been the basis 
of a series of investigations [t| . Since the existence of horizons is probably the most remarkable feature of classical 
gravity that could serve as a link with the quantum description of spacetime, it is important to try and understand 
whether this surface term can arise from the spin-2 field approach. 

C. Plan of the paper 

We will try to address these issues in as straight forward ("dumb") a manner as possible. Section II reviews the 
background material related to spin-2 field and few important results, needed later, are obtained. In particular, we 
introduce a new second rank tensor S ab associated with any matter lagrangian which can be obtained by a well-defined 
procedure. This tensor, in general, is different from the standard energy momentum tensor T ab but coincides with 
the energy momentum tensor for relativistic particle, spin-0 field or spin-1 field. 

In Section III, we start with an action for Einstein gravity for the metric g a \, = rj a b + Xh a b, expand it is in a functional 
Taylor series in h a b and determine the form of the self-coupling at the lowest order. We show that, to the lowest non 
trivial order, the coupling is of the form h ab S a b where S a b is the quantity introduced in Section II. We also exhibit 
the non-analytic nature of the Einstein-Hilbert lagrangian in A and prove that Einstein-Hilbert lagrangian can never 
be obtained by an iteration in A. 

In Section IV we provide a general procedure for coupling the field h a t self consistently to any second rank tensor 
which can be expressed as functional derivative of matter action. This leads to a well defined "rule" for coupling the 
field h a }, to matter fields. We first use it with the tensor S ab we have defined in Section II and show that it leads to 
a generally covariant lagrangian for relativistic particle, spin-0 field or spin-1 field but not in general. We then use 
the same prescription to couple h a b to itself and show that the resulting theory is Einstein's theory (Section IV A). 
Thus, at least in the limited case of spin-2 field interacting with relativistic particle, spin-0 field or spin-1 field [the 
only cases in which we have any observational evidence for gravitational theory!], we have a iterative procedure for 
obtaining the full theory when the self-coupling of h a b is not to the energy momentum tensor. In Section IV B we 
repeat the analysis by coupling h a b using the standard energy momentum tensor T ab . In the case of all matter fields, 
this leads to the standard generally covariant action. But when we use this prescription for coupling h a b to itself, 
we do not get Einstein's theory but a more complicated one which explicitly depends on the background Lorentzian 
metric or the field h a i,. We shall also show (in Section IV C) that previous results, when properly analysed, agree 
with our claims. 

The analysis in this paper goes contrary to the conventional wisdom and Section V discusses the issues which arise 
from this work. 
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II. ACTION AND ENERGY MOMENTUM TENSOR FOR THE SPIN-2 FIELD 



In this section we will collect together the results which are required later. (A more pedagogical description is 
provided in Appendix A; this may be useful since the results are somewhat scattered in the literature |10|\ The action 
for the non-interacting, massless, spin-2 field h ab is built out of scalars which are quadratic in the derivatives d a h bc . 
The most general expression will be the sum of different scalars obtained by contracting pairs of indices in d a h bc dihj k 
in different manner. If we assume that the field equations should be invariant under the gauge transformation: 

hab (x) h ab (x) + d a £ b (x) + d b £ a (x). (2) 

then the resulting expression for the quadratic part of the action can be written in different, equivalent, forms: 

A h = I tfixdahbcdihjk [yf V V fe - ^VV fe + VV V ~ VV V J ] 



~ / d 4 x 



d 4 x 



d t h a a d l h) - d a h bc d a h bc + 2d a h bc d c h ba - 2d a h b b dih l 

\dX a d% - d a h bc d a h bc + 2d a h bc d c h ba 



1 VabK (3) 
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[If we assume that the action is quadratic in the first derivatives and gauge invariant, its form is uniquely given by the 
above equation, except for one very specific four divergence term which can be added. This is discussed in Appendix 
A; see Eq. I|(j5fl .] We shall use the more compact notation: 

Ak = \J d^xdaHAh^M^W™) (4) 

where the tensor M abcl ^ k '(r] mn ) is symmetric in bc,jk and under the triple exchange (a,b,c) <-> (i,j,k) and is given 
by: 

M abci i k (r) mn ) = [r) ai r) bc rf k - r) ai r) bi rf k + 2r] ak rf ,j r) ci - VVy 3 ] symm (5) 

where the subscript "symm" indicates that the expression inside the square bracket should be suitably symmetrized 
in bc,jk and under the exchange (a, 6, c) <-> (i,j,k). In the expression for the action, since 

M abdjk ig mm tiplied by 

dah-bcdihjk, we need not worry about symmetrization and use the expression given inside the square bracket in Eq. JSJ) 
as it is. The gauge invariance of the action leads to the identity 

M abcijk dbd a dih jk = 0. (6) 

To the lowest order, we can couple h ab to other fields by adding an interaction lagrangian of the form (X/2)T ab h ab 
where T ab is some tensor built out of the matter variables and A is a coupling constant. The total action will be: 

Aot = \ f d A x d a h bc d t h ok M ab ^ k {n mn ) + ^Jd 4 x h ab T ab + A matter (7) 

Obviously, only the symmetric part of T ab is relevant for this coupling. The variation of h ab will now lead to the field 
equation M abci i k d a dih ik = XT bc , The condition in Eq.© now implies that d a T ab = 0. Thus the field described by 
the action in Eq. J2J can only be sourced by a conserved, symmetric part of a second rank tensor. The above fact — 
which is, of course, fairly standard — shows the intimate connection between the conservation of the source and the 
gauge invariance of the field. 

It is this conservation law d a T ab — which leads to an inconsistency if we assume that T ab is the standard expression 
for energy momentum tensor for matter fields. When the matter variables are varied, the equation of motion will now 
be affected by h a b because of the (X/2)T ab h ab coupling. But the condition d a T ab = is equivalent to the equations 
of motion for matter field, when it is unaffected by h ab . Hence, in general, it will not be possible to satisfy both 
these conditions and find consistent set of solutions. The conventional wisdom is to attempt to find a consistent 
theory in which the field equations for h ab should imply not the condition d a T ab = but a modified one of the form 
d a (T ab + t ab ) = where t ab is the energy momentum tensor for the spin-2 field. This will require coupling the field 
to its own energy momentum tensor recursively and the hope is to show that — when the recursion is carried out to 
infinite orders — the resulting theory will be Einstein's gravity. This brings us to the question of defining the energy 
momentum tensor for the spin-2 field. 
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For any system described by a Lorentz invariant Lagrangian L( 
with A representing possible tensor indices, one can show that 



i, d a <pA) where (j)A denotes a generic matter field 



db d a (, 



( dL 



\d{d b 4> A ) 



-S b L 



= 



(8) 



when the equation of motion is satisfied 11]. 
tensors of the form 



This allows us to define an infinite number of conserved second rank 



rpba _ 



d a 4>, 



( dL 



\d{d b (f>A) 



- n ab L 



d c tp 



cba 



(9) 



where i/j cba = —^ bca ig an arbitrary third rank tensor anti symmetric in the first two indices, so that d c dbip cba = 
identically. It is possible to choose ip cba in infinite number of ways and still ensure that T ba is symmetric. Thus 
Lorentz invariant field theories possess infinite number of conserved symmetric second rank tensors any of which can 
be legitimately thought of as an energy momentum tensor. For the spin-2 field this prescription gives 



t p<i = l M P b ^ k d l h jk d c 'hb c - rf"L + d c ^ 



(10) 



This non-uniqueness shows that it is not possible to proceed further without making extra assumptions regarding the 
form of T pq . 

One needs to be clear about the different kinds of ambiguity in the definition of T pq . The first ambiguity is in the 
choice of ip cpq . The second ambiguity has to do with the fact that we can add to our Lagrangian a total divergence 
with a undetermined coefficient [as shown in Ea. (|65H of Appendix A]. This changes the form of T pq . Third, the T pq 
defined in Eq. IjlOl) is not gauge invariant. In fact, one can prove a general theorem ^2] that the energy momentum 
tensor for the spin-2 field cannot be made gauge invariant for any choice of ip cpq . This raises serious questions about 
whether the resulting theory after infinite iteration will possess any trace of the original gauge symmetry. If one were 
to be honest, in the sense that no structures other than those sanctioned by Lorentz invariant field theory are to be 
used, then it is not possible to proceed any further and obtain a unique nonlinear theory. 

Let us, however, reduce the standards of honesty and introduce another definition of a second rank symmetric 
tensor (which is based on what we learnt in our general relativity course but we won't mention it!) along the following 
lines: We note that the Lagrangian for any field also depends on the Lorentz metric rj ab , i.e, L = L(<f>A, d a (f>A, Vab)- 
The functional derivative of the action A matter with respect to rj a b will define a symmetric, second rank tensor which 
we can attempt to use. But since rj a b = dia (—1, 1, 1, 1) is a constant {5A/5rj ao ) is mathematically ill defined and the 
functional derivative actually depends on the rule for its definition. We shall see below that several rules are possible 
but let us first consider the conventional wisdom again. 

We begin by noting that, even though we are in flat spacetime, we can use any set of coordinates to describe the 
physics. Let us assume that, in a curvilinear coordinate system we choose, the spacetime metric is jabix). We will 
further assume that the action in the curvilinear coordinates is obtained by replacing i] a b by j a b, ordinary derivatives 
into covariant derivatives and changing the volume element from d A x to d 4 
dependence on j a b which we shall explicitly exhibit by writing it as 

A((j)A,d(j)A,riab) — > Av(4>A, V(/>A,7afc) = / <P 



Thus the action has a kinematic 



,^4>A,lab) 



(11) 



(The subscript V in A-^ is to remind ourselves that, in obtaining this action, ordinary derivatives have been changed 
to covariant derivatives; this will turn out to be important later on.) It is now possible to obtain a second rank 
symmetric tensor T ab by taking the functional derivative of the action with respect to j a b and then setting -f a b = rj a b'- 



SA V = - I d A x V=7 T S^ a b ; T ab (x) 



SA X 



-7 $Jab(x) 



More explicitly, this leads to the energy momentum tensor: 

2 f dLy^f 



T ab (x) 



djab 



\d(d c j ab ) 



(12) 



(13) 



7=17 

The procedure described above provides one possible prescription for obtaining T ab in flat spacetime. Note that jab 
for us is purely a bookkeeping device and, in the end of the calculations, we shall set j a b = r\ a b- 
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It is rather surprising that this definition of T ab is routinely used in field theoretic approaches to gravity (like, for 
example, in 0) as though it has nothing to do with curved spacetime. This attitude is incorrect. The variation of j a b 
to 7 a fc + Sj a b for arbitrary choices of Sj a b takes one from flat spacetime in curvilinear coordinates to genuine curved 
spacetimes. (Just varying the coordinates in flat spacetime will only have 4 function degrees of freedom while we need 
10). The evaluation of the functional derivative in Eq. requires the strong assumption that the action in Eq. Hll|) 
is valid in arbitrary curved spacetime with metric jab- We have come a long way from the basic concepts of Lorentz 
invariant quantum field theory in flat spacetime. 

To use this definition for the spin-2 field, we first write the action in Eq. in arbitrary curvilinear coordinates 
with metric j a b, using our rule in Eq. i|ll|) as: 



A 



^ f / 7" T7 u T7 h r ^Ai^fic^jk ni ai ni bj n ck , n^,,ak^,bj „ci r)^,ak^bc„ij~\ 

= -ax v^7 vahbcVihjk [1 7 7 — 7 7 7 + l l 7 7 ~ l l 7 7 J 



(14) 

In this expression, the covariant derivative operator V is defined with respect to the metric j a b and involves the first 
derivatives d a jbc- Varying this action with respect to -f a b will lead to a symmetric energy momentum tensor for the 
spin-2 field when we use the prescription in Eq. Q13[). The actual expression for this tensor is fairly complicated but — 
fortunately — we do not need it. We, however, stress the following fact: Since [dL/d(d c j a b)] involves first derivatives 
of h a bt the tensor T ab will involve second derivatives of h a b. It can again be shown by detailed algebra that this 
T ab is indeed of the form in Eq. 11U|) for a specific choice of ip cpq ; thus, our rule chooses one out of many choices in 
Eq. (|10l) . We can now attempt to obtain the non linear theory by coupling h a b iteratively to the energy momentum 
tensor defined by Eq. (| 1 31) - We shall show in Section IV B that contrary to popular belief, the resulting theory is not 
Einstein's theory. 

There is, however, a more important issue which needs to be raised as regards the procedure used to obtain Eq. 1)13(1. 
What we have done is essentially to introduce the curvilinear metric j a b into the matter action (which was originally 
defined in flat spacetime Cartesian coordinates) by a particular rule and then evaluate the functional derivative in 
Eq. ifTBT) . At the end of the calculation, we set j a b — > r/ a b- The rule in Ea. (|TT|) is strongly motivated by general 
covariance and, of course, leads to a generally covariant matter action in the curvilinear coordinates. But since we do 
not have the right to assume general covariance (and only Lorentz invariance), the rule we have specified is only one 
among many possible ways of introducing j a b into the matter action. 

To bring this point sharply into focus and to derive some important consequences in the coming Sections, we shall 
introduce another rule which leads to the definition of another symmetric second rank tensor. To do this, we will 
construct a modified action in the curvilinear coordinates by replacing r\ a b by j a b and changing the volume element 
from d x to d x^J — 7 but without changing ordinary derivatives into covariant derivatives. The action again acquires 
a kinematic dependence on j a b which we shall explicitly exhibit by writing it as 



A((f> A , d(j>A, Vab) -> A d ((j>A, d(j) A , lab) 



-7 L d {4> A ,d(j)A,lab) 



(15) 



(The subscript d in Aq is to remind ourselves that, in obtaining this action, ordinary derivatives are retained as they 
were.) We can again obtain a second rank symmetric tensor S a b by taking the functional derivative of the action with 
respect to jab and then setting j a b = r) a b- 



SA a = 



ab { \ 



SAf 



"7 hab{x) 



(16) 



More explicitly 



S ab (x) 



d^ab 



= 2 






- 1=V 


dlab 


7=>) 



(17) 



This procedure provides another possible prescription for obtaining S ab in flat spacetime. Once again, the j a b is 
purely a bookkeeping device and, in the end of the calculations, we shall set j a b — r) a b- 

Both S ab and T ab are Lorentz invariant tensors. In general, there are two crucial differences between these two 
tensors T ab and S ab : (i) The T ab is obtained from a generally covariant lagrangian and hence is generally covariant; 
the S ab need not be generally covariant (ii) The T ab satisfies the identity V a T ah = since it arises from an action 
which is a generally covariant scalar; the S ab need not satisfy this identity. Having said these, one must note that 
these two tensors are identical whenever the action does not depend on the derivatives of the metric. For spin-0 field, 
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spin-1 field and for a relativistic particle, the generally covariant action in Ea jlljl is independent of the derivatives 
dalbc of the metric tensor. Hence, in all these three — physically important cases — the two definitions lead to 
identical energy momentum tensors: T ab — S ab . Even our apparently non-covariant definition will lead to a generally 
covariant energy momentum tensor which satisfies the condition W a S ab = 0. 

For the spin-2 field these two definitions differ. To find S ah , we need to start with the action in the form 

jd a h bB d i h jk M abciS \'y mn ) 

jd a h bc dih jk [ 7 <« 7 bc y fc - y» 7 V fc + 2 7 a V J 7 " - 2 7 afc 7 b V J ] (18) 

which differs from that in Eq.JHJ by the fact that we have not changed cVs to V a 's. The tensor S ab can now be 
calculated using Eq. lfT7|) : 

d a h bc dihjk (19) 

Again it is possible to write down the explicit expression for this but, fortunately, we will not need it. But it is obvious 
from this expression that this tensor is quadratic in d a hb c and does not involve second derivatives of h ab . We shall 
see in the next section that it is this object which governs, through a term S ab h ab in the Lagrangian, the coupling of 
gravity to itself at the lowest nontrivial order. 




S pq (x) 



2 dL^P 

^1 djab 



Sy^yM" 3 fc ( 7 mn ) 



<9 7 
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III. SNEAK PREVIEW: REVERSE ENGINEERING OF EINSTEIN'S THEORY 

We want to obtain Einstein's theory by starting from the action for a spin-2 field in flat spacetime and coupling it 
to some kind of energy momentum tensor for h ab and iterating the process. Given the ambiguities in the definition of 
the energy momentum tensor described in the last section, it makes sense to do it the other way round and identify 
the correct form of the tensor to which h ab couples to in Einstein's theory. This exercise is straight forward: (i) Start 
with an action functional A g [g ab ] which leads to Einstein's field equations for the metric tensor g ab . (ii) Define the 
spin-2 field through g ab — r) ab + Xh ab , where h ab has the dimension (length)" 1 and A has the dimensions of length, 
(iii) Substitute g ab = i] ab + Xh ab in [5at] and expand in a Laurent-Taylor series in A [or, which is the same thing, 
do a functional Taylor series in h a b\. Now, if our ideas are correct, two things must happen: (a) the lowest order term 
should give the action functional for spin-2 field in flat spacetime with a suitable choice for A; (b) the next order term 
will have a lagrangian of the form Xh ab K ab and we should be able to read off K ab . We will carry out this exercise 
and then comment on various issues. 

The conventional action principle for general relativity is the Einstein-Hilbert action given by (with A 2 = 4irG) 

A EH = 7^2 / RV^gd^x ee J ^[V^iquad - 9jP j ] ee A quad + A sur (20) 

where 

/"'(i;;,,!;, KA) (21) 



L, 



quad 



and 



p c = V=g (g ck rTm - 9 ik r c ik ) = V=g (g ac g Jl - g ta g cj )d igac (22) 



The quantity L qua d is what is usually called the T 2 lagrangian and is quadratic in the first derivatives of the metric. 
The term diP 1 integrates to a surface term and is usually ignored while deriving the field equations by assuming that 
"suitable" boundary conditions can be imposed. A more formal route is to add a suitable boundary term to cancel 
this, thereby essentially reducing the action to one based on L qua d- 

Classically, there is no way of deciding whether Aeh or A qua d is the "correct" action, since both lead to the same 
field equations. Let us first consider A qua( i and determine to what second rank tensor it self-couples at the lowest 
order: Since T ~ dg, one can express A quac i as a quadratic expression in d a gbc- Straightforward algebra gives this to 
be of the form: 

Aquad = / d*x V=9 d a g bc d i9jk [g ai g b W k ~ g ai g bj g ck + 2g a k g bj g ci - 2g ak g bc g i3 } 

= 7^7 J d i x^d a g bc d i g jk M abci ' k (g mn ) (23) 
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with the same functional form as the M abc ^ k defined in Eq. J3J ! This is a miracle and all the results of this paper 
are essentially an exploitation of this miracle. This result means that, if, after obtaining the flat spacetime action for 
spin-2 in Eq.Q, we have (i) "just replaced" r] a b by g a b in M abc ^ k and d A x by d A x^—g and (ii) used g a b = Vab + Xh a b 
in the derivatives d a hb c , we would have got the correct Einstein's theory. All the (infinite) iterations are required only 
to understand why this is legal. 

Let us now proceed with the original programme of reverse engineering this action to find out what it couples to at the 
lowest order. This is quite straightforward. Substituting g a b = r/ a b+Xh a b in Eq. I]23p. we get d a gbcdigjk = X 2 d a hb c dihjk 

abcijk 



exactly. The Taylor series expansion of \J—gM 



gives: 



yf^g~M ahci i k = M ahci i k (r\ mn ) + X 

L d 9 Pq 
Putting them together, we get the expansion: 

Aquad = j f d A xd a h bc d l h ]k M ab ^ k {ri mn ) + ^ { d l xd a h bc dih jk 



f abcijk ^gTnn 



h PQ + 0(\ 2 



r abcijk ( n mn 



{g mn ) 



ipq 



o(x 2 



(24) 



(25) 



But the integrand of the second term contains precisely the quantity we defined in Eq. Ijl9|l . Using it, we get the final 
answer: 



-4, 



quad 



d A xd a h hc d 1 ,h Jk M ah ^ k {T ] mn ) 



x 



d 4 xS pq h 



i'<i 



0{X 2 ) 



(26) 



with S pq given by: 



gpq 



jabcijk ^mn 



dlpq 



dJibAh 



jk 



(27) 



7 = 7) 



We, therefore, have proved that the coupling of gravity to itself, at least to the lowest order, is to a strange beast, 
defined in a noncovariant way. To 0(A), the field h a b couples to a quantity which is quadratic in the first derivatives, 
dahbc: of the field and does not couple to an object which has second derivatives of the field. The standard energy 
momentum tensor T pq defined in Eq. (|13|l involves second derivatives of the field and hence a naive coupling of the 
form T pq h pq will not match with what we have found by explicit computation above. (The proof in Eq. i|26|) . of 
course, gives the result only to the lowest order and one needs to know whether gravity consistently couples to S ab 
[as defined by Eq. 119|) ] to all orders in the coupling constant. It should be intuitively obvious that it will, but we 
shall provide an explicit proof in Section IV A.) 

It is, anyway, easy to understand that the self-coupling term in gravitational lagrangian cannot be of the form 
\h ab Tg b where Tg b is the energy momentum tensor of the graviton. The reason is the following: A term in the 
lagrangian proportional to Xh a bT ab where T ab is due to external matter fields (assumed to be independent of h a b 
to this order), will lead to the equation of motion of the type d 2 h = XT. Hence a coupling of the type Xh a bT ab is 
equivalent to requiring the source being T a b- Consider now the coupling of gravity to itself through a term of the type 
Xh a bt ab (h) where t a b is some tensor which explicitly depends on the graviton field h a b- When this term is varied with 
respect to h a b to get the equations of motion we will obtain two terms: t ab + (dtij/dh a b)h 1 ^ both of which will act 
as a source to gravity at next order. If we want the source to be the energy momentum tensor of graviton field, TR, 
then the coupling cannot be of the form h a bT ab (h)c since this will lead to the wrong source. What we find is that the 
coupling in the lagrangian should be to S ab if the source of the gravity is to be the energy momentum tensor to the 
lowest order. (In this paper, we shall use the terminology, " A is couples to be B , \ if the lagrangian at the relevant 
context has the term AB. This does not necessarily mean that B acts as a source term in Euler Lagrange equations 
when A is varied, since — in general — B could be a functional of A.) 

So far we have obtained the lowest order self-coupling from the L 2 action and we will now turn our attention to 
the surface term. We will now prove a strong result: It is impossible to obtain the Einstein- Hilbert action, especially 
the A sur term, by starting from an action for spin-2 field which is quadratic in the first derivatives and iterating in 
powers of X. To see this qualitatively, note the structure of the Taylor series expansion for A sur ; symbolically: 



-4., 



1 

A 2 



d 4 x d[Xdh{n + Xh + X 2 h 2 + 0(A 3 ))] 



d 4 xd[ndh}+ / d i xd[hdh] + X / d A x d[h 2 dh] + 0(X 2 ) (28) 



(The Lagrangian in A sul - is linear in first derivative of gab- 
type term which — on multiplication by the pre-factor A~ 2 



The substitution g a b — rj a b + Xh a b will lead to a Xdh 
— will give rise to the lowest order term which scales as 
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A 1 ). The third term of the above expansion, which is 0(A) and has the structure d(h 2 dh), can be combined with 
the coupling term (X/2)S pq h pq in Eq. which is also 0(A). If we write: 



d{h 2 dh) ~ h[hd 2 h + {dhf 



(29) 



it might seem that the S pq in the coupling term (X/2)S pq h pq in Eq. Il2(ill changes with a contribution from the second 
derivatives of the field. One might wonder whether this will help us to get A sur by using the coupling (X/2)T pq h pq 
and exploiting the second derivatives d 2 h in T pq . This idea, however, will not work. The term we need [0(A) term] 
is the third term in the Taylor series expansion in Eq. l|28[l and to get Aeh as the final answer, we must obtain first 
two terms as well. The two leading terms in Eq. I|28|l are of 0(1/A) and 0(1) and hence cannot be obtained by the 
iterative process by coupling in ascending powers of A. 

More explicitly, the two leading terms in A sur in the Einstein-Hilbert Lagrangian are: 



^■Jd^dadblh^-^h^ 



4A 



d A x d r 



hid c h - hid k h kc - h kc d k h - h ab d c h ab + h hc d k h\ + h ak d k h c a 



0(A) 



Thus the leading term in Einstein-Hilbert Lagrangian (the term in the first line of the above equation) 



A. 



~J d 4 xd a d b [h ab - V ab h^ 



(30) 



(31) 



is non analytic in A when expanded in g ab = rj a b + Xh ab . Note that one cannot get out of this non analyticity by 
cheap tricks like rescaling h ab to Xh ab . The dimension of the genuine graviton field and the form of the zeroth order 
lagrangian uniquely fixes the scaling of h ab and there is no freedom for dimensionful scaling left. (In fact, it turns out 
that this non-analyticity is vital for the interpretation of the surface term as horizon entropy in semiclassical gravity; 
so it is not a trivial issue. This is discussed in detail in a separate publication [13J ) For us, the importance of the 
non-analyticity lies in the following fact: If one starts with the quadratic spin-2 graviton action and iterate with self 
coupling, it is impossible to obtain R^J—g, since it requires obtaining a piece non-analytic in A. One may wonder how 
previous "derivations" obtained this piece. This was added by hand; since all the terms in A sur are four divergences, 
any part of it can be added by hand without affecting the equations of motion and this is what was done. (We will 
discuss this in detail, in the context of the derivation by Deser [j| in Section IV C.) 



The Lagrangian in the leading order surface term in Eq. (|31[1 . viz., L s 



d a d b [h — rf h\], which has some 



interesting properties. First, it is gauge invariant under the transformation Eq.Q, as can be easily checked. (This, 
of course, means that one cannot set this term zero by a gauge choice). Second, it provides a simple counter example 
to the belief that if a functional is invariant under infinitesimal gauge transformations [that is, under Eq.JSJ), with £ 
treated as first order infinitesimal], then the expression will be generally covariant. This belief originates from the 
fact metric tensor transforms like in Eq.JSJ under infinitesimal coordinate transformations (x l — ► x l +£ t t (x)) and one 
thinks of [erroneously] the finite coordinate transformations as arising from "exponentiating" the infinitesimal ones. 
Explicitly, the functional 



F 



-gd l9jk {g jk g li -g lk g ji )] 



(32) 



is clearly not invariant under arbitrary coordinate transformations. But if we take g ab = J] a b + Xh ab , then the expression 
for F becomes, to linear order in h, 



F Vm » Xd t [d t h jk (rf k r, H - rj'V)] = -XdA 



rf l h] 



(33) 



which is the same as the integrand in Eq. (|31|l . This expression, however, is gauge invariant under the transformations 
in Eq. . This shows that it is possible to have scalars which are invariant under infinitesimal gauge transformations 
but are not generally covariant. 



IV. FROM GRAVITONS TO GRAVITY: GENERAL PROCEDURE 



After the sneak preview of the results to come, we shall now return to the original task of coupling the spin-2 field 
to matter, as well as to itself, self consistently to all orders. We shall start with the issue of coupling the spin-2 field 
to other matter fields self-consistently to all orders to see how an externally specified, h ab {x) affects the dynamics of 
(j) A {x). 
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Consider a field (j>A(x a ) described by a Lagrangian density L((j>jH,d(j>A,r]ab) in flat spacetime, in the Cartesian 
coordinates in which the metric is r] a t = dia (— 1, 1, 1, 1). The index A formally denotes all the indices the field carries 
depending on its spin; we will assume the field is bosonic for simplicity. To couple h a b to this matter field, we need to 
first find a suitable, second rank tensor field K ab , defined in terms of the matter variables. In Section II, we introduced 
two such tensors (among infinite number of possibilities) in Eq. 1|13|) and Eq. 117|) . Both can be generically expressed 
as functional derivatives in the form: 



SA = ^J d i x^/K ab 5 lab - K a \x) ee 



5A 



-1 S"fab(x) 



(34) 



If we use Aq = A-^ijpA, V^a, "fab) (with covariant derivatives), we get the conventional energy momentum tensor 
K ah = T ab . If, instead, we use A = Ai)((j)A,d(j)A,Jab) we get K ab — S ab which is a hybrid object, except for the 
relativistic particle, spin-0 field or spin-1 field for which the definitions coincide and T ab — S ab . For most part of our 
algebra below, we need not specify whether we are using A-y or Aq and we will use the generic symbols A and K ab 
to stand for either definition. 

We now want to couple the second rank symmetric tensor field h a b to K ab . To the lowest order, this is done by 
changing the action from Aq to A<i = Aq + A\ where A\ is chosen such that: 



8A X 



A 



<f xiV-^ K^ix^SKbixx) = A / d x 



SA 



Sjab(xi) 



5h ab (x 1 ) 



(35) 



1=V 



where A is a coupling constant. To the lowest order, K ab is independent of h a \, and we can integrate this to obtain 
the action: 



A< X = Aq + Ax = Aq + A / d 4 xx 



sa 



.Sjab(xi) 



7=') 



h ab (xi) = A + A / d 4 x ls /^K a, '(x ) CU.) 



The addition of this coupling will, however, change the definition of K ab , since the second term A\ contributes to 
K ab via Ea. (|34l) . To take this into account, we need to add a term A-i in a manner similar to what we did in Ea. (|35|l : 
that is, we need to impose: 



SA 2 A / d A x 2 
Using the form for A\ in Ea. (|3l)|) we get 



SA 1 



SA 2 



d x\d x 2 



hcd{x 2 ) 



5 2 Aq 



Sh cd (x 2 ) 



7=17 



_57cd(x 2 )S"fab(xi) 

Integrating, we get the second order correction to be 



h ab { Xl )5h cd {x 2 ) 



A 2 

A 2 = — I d 4 xid 4 x 2 



S 2 A 



S"fcd(x 2 )6"f ab (xi) 



h ab {xi)h cd {x 2 ) 



(37) 



(38) 



(39) 



7=1 



It is obvious that this term will bring about another correction etc. The sum of the infinite series of terms in the 
action will be 



-4. 



00 \ n r 



. d X r , 



S n A 



hab{x\) hij(x n ) 



(40) 



which is just a functional Taylor series expansion leading to: 



A x = A (j ab + \hab) L 



(41) 



An alternative way of obtaining the same result is note that, at every order, we have the recurrence relation, similar 
to Eq. and Eq. (J3SJ: 



SA n+1 = A 



J d i x 1 



SAn 
5-fab{x 1 )_ 



ShabiXx) 



(42) 



7=7) 
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which is the same as: 



8A n+ i 
\5h ab (x) 



6A n 



Summing both sides to to infinite orders, we find that Aoq satisfies the relation 



SA C . 



\Sh ab (x) 



SA r _ 



(43) 



(44) 



which has the general solution given by Eq.l^TJ- This is our key result. 

Thus we can consistently couple a field h a b to K ab by the rule given in Eq. (|41|l : This allows us to compute the 
effect of an external h a b on the system if we insist that the external field consistently couples to a tensor K ab which 
can be expressed as in Eq. I|34(l . Since the curvilinear metric 7a b was introduced only as bookkeeping device to allow 
for variation of the action, the final action is given by replacing j a b by rj a b at the end of the calculation. Note that 
there is subtle difference between Eq. <|41|) and the expression obtained by replacing n a b by (rj a b + Xh a b) in the original 
action. The latter one will miss, for example, the -^/det^l — > ^/ det 1 77 + h\ kind of factors. We need to introduce 7a & 
in order to provide a placeholder in the final expression. 

Here comes the parting of ways. If we had chosen A — A^j{4>a,^ 4>A,lab) (with covariant derivatives), then 
K ab = T ab is the standard energy momentum tensor and our result gives the final matter action to be: 



Aao = A ((j} A ,y(g)4'A,gab) 



(45) 



where we have used the abbreviation g a b = n a b + Xh a b and V is defined with respect to this metric. This is a generally 
covariant matter action in the spacctime with metric g a b and agrees with all the text book results. It should, however, 
be stressed that this cannot be considered a derivation of general covariance of matter action when self-consistently 
coupled to spin-2 field. This is because we made a rule for finding T ab which has general covariance with respect to 
curved spacetime built in as an assumption. All that has been shown is that this extends to an interpretation of g ab 
as a metric tensor and only the combination g a b = rj a b + Xh a b is relevant for matter sector. 

Since the final matter Lagrangian we obtained is generally covariant, its energy momentum tensor has zero covariant 
divergence, leading to: 



5A C 







where V is the covariant derivative operator corresponding to g a b = r) a b + Xh a b- 

On the other hand, if we had chosen A = Aq(4>a, d4>A,Jab) (without covariant derivatives), then K ab 
beast we have introduced in Section II D and our result gives the final matter action to be: 

Aoo = A (4>a, dip a, 9ab) 



(46) 



S ab is the 



(47) 



In general, this is not a a generally covariant matter Lagrangian since we have not replaced partial derivatives by 
covariant derivatives. The metric appears only through \J—g factor and by the replacement of r\ a b by g a b- This, of 
course, does not matter for the lagrangians of relativistic particle, spin-0 field or spin-1 field. In all these cases this 
prescription does lead to a generally covariant lagrangian though this is not by design. The Ea. (|46[) also holds for 
the spin-0 or spin-1 field or for the relativistic particles but not in general. (We shall comment on this in Section V.) 



A. How to obtain Einstein gravity from the spin-2 field 

To see where all this is leading to, let us consider next the real issue: that of coupling the graviton field to itself. In 
our approach this is ridiculously simple. We merely use the fact that the analysis leading to Eq. I|41|) was completely 
independent of the form of Aq as well as the nature of the fields 4> A . Hence we can use the same prescription when 
if> a is the second rank symmetric tensor field h a b itself. 

The key ambiguity, of course, is whether we want to use A\? and couple to T ab or use Aq and couple to S ab . Let 
us use the second procedure first: The Aq for the graviton field is obtained from the first line of Eq. Q by replacing 
rf b by 7 afc in M abci i k {rj) and multiplying by \/ = 7 . This is given by Eci. (fT£|) : 

Aq = - A J d i X^d a h bc d i h jk M ab ^ k ( 7 mn ) 

= \J d A x^ld a hb c d t h 3k [ 7 <» 7 V fc - y» 7 V* + 2 7 a V J 7 " " 2 7 a W J ] (48) 
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Our prescription now requires that Aoc for the field h ab is obtained by replacing ^ ab by g ab = Vab + Xh ab . This leads 
to the action 

Ao ° = i J d A x^d a g bc d l9jk M ab ^ k (9 mn ) 

= j "d 4 xV=g d a g bc d i9jk [g al g bc g jk - g az g bl g ck + 2g ak g b ^g cl - 2g ak g bc g^] (49) 

which is precisely the T 2 action in general relativity. The variation of this action will lead to Einstein's field equations 
in vacuum and we have achieved our goal. 

As an aside, we note that the above approach will work with any generic action functional for spin-2 field. A 
different choice of the spin-2 action functional will lead to a different final theory. This aspect is briefly discussed in 
Appendix B since it is irrelevant to the main issues of this paper. 



B. The failure of the conventional procedure to lead to Einstein's theory 

Let us see what happens if we follow the conventional procedure. For this, we need to start with the Asj of graviton 
action and couple to the standard T ab . This action is given by Eq. I)14|l: 

A V = i j d^X y/^VahbcVihjk M ab ^ k (j mn ) 

= \ J d^x^Vah^^ [ 7 « 7 V fe - 7 <W* + 2 7 a - 2 7 afe 7 bc 7 ^ (50) 

Our prescription now requires that is obtained by replacing j ab by g ab = i] ab + Xh ab . This leads to the action 
Aoo = \ J d^x^—gV a h b ^ l h, k M ab ^ k {g mn ) 

= \ J d^x^g-VahcV^k [g al g bc gi k - g al g b >g ck + 2g ak g b ig™ - 2g ak g bc g^} (51) 
where the V operator is with respect to the metric g ab . Now, since Vigjk = 0, we get: 

V a h bc = jV a [g bc - T] bc ] = -j^aVbc = -jFlaVic + KaVib] 

= -\[g P ^ P ba(g lc - AM + (b «-> c)] = -j[T pba (SP - Xg pl h lc ) + (b^ c)] (52) 

and a similar expression for Vihj k - Since these are multiplied by M abci: > k (g' mn ) which is symmetric in (b, c), we 
can ignore the (b <-> c) term etc. We thus obtain 

VahcVihjk M abcijk (g mn ) = ^[T pba T qkl {5 p c - Xg pi h lc )(S] ~ Xg l ' 1 h lj )]M ab ^ k (g mn ) (53) 

Of the four terms which arise on expanding out the product, (S P — Xg pl hi C )(Sj — Xg lq hij), the first term can be 
transformed, again using the symmetry of M abci: > k (g mn ) in (b, c), (i, j) to give: 

^T cba T jki M ab ^ k (g mn ) = ^M ab ^ k (g mn )[T cba + T bca }[T jki + T kji ] = ^M ab ^ k (g mn )d a g bc d t g jk (54) 

This is precisely the lagrangian term in Einstein's theory, in the form of T 2 action [compare with Eq. 149J1 ]. Unfor- 
tunately, there are three more terms in Eq. (|53[) of the TThh and TTh form. They do not vanish, they are not total 
divergences and they depend explicitly on h ab . Thus the action functional obtained by coupling the spin-2 field to 
the standard energy momentum tensor is not that of Einstein's theory. 

The manner in which we led the reader to the result, it should not come as a surprise. We have, in fact, shown in 
Section III explicitly that the standard action for Einstein's theory does couple to S ab in the lowest order. So clearly, 
we will not get Einstein's theory if we force feed the conventional energy momentum tensor T ab on to the spin-2 field. 
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C. Comments on the previous work 



Finally, we shall discuss how the previous "derivations" escaped this problem. As described in section I, none of 
the previous derivations other than that of Deser [j| actually performs any iteration. All the rest of them tacitly or 
explicitly brings in general covariance for the gravity sector, after which it is trivial to obtain R^—g as the Lagrangian; 
hence, we need not discuss them any further. Deser does perform the iteration using the first order form of gravity 
and using \J~gg ab as the chosen variable. Let us first summarize this approach briefly in a slightly different language 
to bring the essential ingredients in to focus. 

We first note that the standard Einstein-Hilbert action can be expressed in the first order Palatini form using the 
variables T bc and f ab = \f—gg ab as: 

AEH = i / d V ab i?a 6 (r); f ab ee ^—gg ab (55) 

Varying f ab and T bc independently in this action will lead to standard Einstein's equations. If we substitute f ab = 
rf + Xq ab into this action (without any approximations), then the lagrangian becomes: 



^EH 



L [(v ab + M ab )K b F) + R%™ (56) 



4A 2 ' 

where R^ b (T) and R® b (T) are the linear and quadratic parts of the Ricci tensor 

R^ b (T)=d c T c ab -d b T a ; R^ b (T)=T c r c ab -T c ad Tt c ; T c = P ci (57) 
Expanding out the product in Ea. (|56() . we get four terms which we will group as: 

L EH = L +L 1 + L 2 = J_tf*i&(r) + J^i&tT) + Xq ab R^ b (T)) + ^Xq ab R^ b (T) (58) 

We now notice something remarkable. 

• The first term i] ab R^ b = di(T l — r JJ ■) is a total divergence. Let us assume we are allowed to drop this term. 

• The second term has the piece [r/ ab 'R^ b (T) + Xq ab R^ b (T)] which is essentially equivalent to the zeroth order 
action for spin-2 graviton (in the second order formalism) plus a very specific four divergence term. [This will 
be apparent if we write the term qR L ~ qdT ~ d(qT) — Tdq; we find that it has a quadratic term plus a very 
specific total divergence term d(qT)]. Let us assume we are allowed to start with this very specific term as the 
lowest order graviton lagrangian. 

• Granted these two wishes, we see that one can obtain from this L\ term, the tensor t ab = (5Li/Srj ab ) — 
{l/AX 2 )R%{T) , purely formally. If we think of this as the energy momentum tensor for the graviton, then 
the next order coupling should be Xq ab t ab = (l/4A 2 )[Ag afc i?^ b (F)] which is precisely the last piece L 2 in the 
Einstein-Hilbert lagrangian ! What is more, this last term is independent of rf b so (8L2/ 8r) ab ) — and no 
further iterations are required. [Since we are working with \f-~gg ab rather than g ab , the variation actually gives 
Tat — (l/2)3afcT' rather than the energy momentum tensor itself; but this is irrelevant to our discussion.] 

When one attempts to do this properly, one faces an important issue: As we have said before, one cannot really 
define things like (b~L/8r) ab ) where ■q ab is the Minkowski metric. So we first need to write L\ in curved spacetime 
with a metric j ab and compute variations with respect to this metric. What do we do to the d a in the definition of 
Rab(T) = dc^ab ~ dbT a l Let us suppose, we change them to V Q for the metric j ab . Then the t a b computed from the 
functional derivative will pick up additional terms. The action L\ in curved spacetime is 

£i(7) = 4^[V=77 ob ^ b (r) + Xq ab R L ab {Y)\- r b = V^l ab (59) 

where the R ab is now evaluated with partial derivatives replaced by covariant derivatives with respect to "f ab etc. The 
variation of this Lagrangian with respect to j ab gives (see eq. (8) of 0) 



tab ~ 4A 2 



RW) + X °ab (60) 
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where the second term a ab arises from the variation of the da'ybc terms, because we have changed d a to V a . Its explicit 
form 



<?ab = C 



^lab [ qfK 3 - \q T c ) - 2q l c T [ab)l - 2q l {a T clb) + 2q l [a T b)cl - q ab T c + 2q c(a r b) 



(61) 



is irrelevant to us. We are now in trouble since the next order coupling t ab q ab will now have an unwanted term q ab a a b, 
in addition to the term we want (proportional to q ab R^ b )- Deser simply drops this term saying (see his comment after 
eq. 9 of 0): "Note that we have not added the full Wt^, but rather used the simple part of t^ v only" without any 
additional justification! Then, of course, one gets the Li as the next term and iteration stops there. 
It should now be obvious that Deser's derivation requires the following implicit assumptions: 

• One should drop the L$ term in Eq. I|58|l unceremoniously, saying it is a total divergence [Deser does this in 
going from eq. (2) to eq. (4) in |5(]. It is precisely this term, which, in second order formalism, has the non 
analytic behaviour (1/A) and is displayed as the first term in Eq. (|30[1 . Sure, it is a four divergence, but one 
can never get it from graviton's quadratic action and one needs to add and subtract this term, at will, to get 
Einstein-Hilbert action. 

• One should start with L\ which is not the graviton action that is quadratic in the derivatives of the field (that a 
particle physicist would have written down from first principles) but the one with a very specific total divergence 
added to it. This is precisely the 0(1) in Eq. (|3U|) . There is no way anyone could have guessed this specific 
total divergence term without knowing the final answer! 

• One should drop the terms in t ab which arise from varying d a ^/ bc . But this is precisely the same as using our 
quantity S ab ! Or, rather, not changing d to V when one takes the graviton action from flat spacetime to curved 
spacetime. So in real terms the two derivations match mathematically and our conclusion stands: Gravity 
self-couples to S ab ; not to T ab . It just was not realized before. 

If we are not attempting to derive Einstein's theory from the spin-2 theory but only want to reinterpret it in the 
field theoretical language, then one may be willing to live with the first two issues mentioned above. One can use 
the hindsight gained from general relativity and add and subtract four divergences at will to the action, to obtain 
the necessary final form. (We must then admit that the Venusian physicists whom Feynman keeps referring to in £| 
would never have got there.) But the third issue is not a matter of opinion or point of view; what quantity gravity 
couples to in becoming nonlinear is a well defined mathematical question. Our analysis — and proper interpretation 
of previous work — gives an unconventional answer. 



V. CONCLUSIONS 



We have shown that it is not possible to obtain the Einstein-Hilbert action starting from the standard graviton 
action and iterating in the coupling constant. This is because of the existence of the total divergence term in the 
Einstein-Hilbert action which is non-analytic in the coupling constant, when expanded in terms of the graviton field. 
This result is crucial because a series of previous investigations [9| have shown that the surface term is vital in the 
thermodynamics of horizons and in semiclassical gravity. In fact, I started this investigation to understand how the 
surface term — and hence, possibly the entropy of horizons — can be interpreted in terms of graviton field in a 
Minkowski background. The result shows that one simply cannot understand the surface term in a standard field 
theoretical language, u sing the graviton field. There is more to gravity than gravitons and this will be elaborated in 
a separate publication [l3| . 

In a strictly classical theory, what matters is the equation of motion and not the form of the action principle. Hence, 
the fact that we can not get the surface term in Einstein-Hilbert action is not of concern if we are only interested in 
the Einstein's equations. Our analysis shows that it is indeed possible to obtain the quadratic T 2 action (and thus 
the Einstein's equations) by starting from the the graviton action and iterating on the coupling constant. But to do 
this, we need to couple h ab to a second rank tensor S ab which is different from from the standard energy-momentum 
tensor Tq of the graviton. Indeed, as we explained in Section III (see the discussion after Ea. (|27|l 'l. if the source of 
gravity at each order of iteration has to be the energy-momentum tensor of the graviton evaluated at the previous 
order, then the coupling in the Lagrangian cannot be of the form h ab TQ b since the h ab dependence of the Tq 1 will lead 
to an extra term on variation. A term in the lagrangian of the form h ab S ab does lead to the energy-momentum tensor 
as the source of gravity. Identifying the nature of S ab and bringing it into focus has been one of the results of this 
paper. 
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If we were only interested in pure gravity, this would have been the whole story. But, in that case, it is an 
unnecessary exercise. The linear spin-2 field, uncoupled to anything, is a perfectly consistent theory and we need not 
try to couple it to itself. So the whole exercise has meaning only when we have both matter and spin-2 field and we 
try to couple them consistently. Then we need to assume that the spin-2 field couples to itself through S ab while it 
couples to matter through T ab . This assumption will lead consistently to Einstein's theory and seems to be the most 
viable option, if we want to obtain standard gravity coupled to matter, starting from the graviton action. [Of course, 
in a world made of a spin-2 field coupled to matter made of only relativistic particles, spin-0 fields and spin-1 fields, 
one can assume that all the coupling is through S ab ; this is because for matter made of these constituents, S ab = T ab ]. 

Two facts need to be borne in mind as regards this option. First, we do not know anything about the coupling of 
spin-2 field to itself except through standard gravity; and the analysis in Section III shows that gravity does couple to 
itself through a term h a bS ab . Second, there is no conflict with principle of equivalence even though the self-coupling 
term is h a bS ab while the coupling to external source is through h a bT ab . What matters for principle of equivalence is 
the fact that the source for gravity is always the energy-momentum tensor This is indeed assured in our approach 
and — as have been stressed several times — this requires a self-coupling term of the form h a bS ab in the lagrangian. 
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Appendix A: Physics of the spin-2 field - brief review 



In this Appendix we will briefly review the theory of the spin-2 field and collect together different results which 
are required later. (This is done especially since I could not find a a convenient source for pedagogical discussion of 
spin-2 field; the results are somewhat scattered in the literature (lfH). 



A. Action functional for the spin-2 field 

The action for the non-interacting, massless, spin-2 field h a b is built out of scalars which are quadratic in the 
derivatives d a hb c - The most general expression will be the sum of different scalars obtained by contracting pairs of 
indices in d a hb c dihjk in different manner. Since this product is symmetric in (b,c) and (j, k) and also under the 
interchange (a, b, c) — > (i,j,k), it is easy to figure out that, a priori, seven different contractions are possible. For 
example, if a is contracted with i, then there are two possibilities for contracting b (with either c or with j; contracting 
b with k is the same as contracting b with j). These contractions will lead to the terms cid a hb c dihjk'n a ' l "r] bc ri^ k + 
c-id a hb c dihj} z , rf l, n b ^ -rf k — Cid a h b d a h^ + C2d a hb c d a h bc in the lagrangian with as yet undetermined constants [c\,C2). 
For brevity, we will denote these two terms symbolically as (ai, bc,jk), (ai, bj, ck). Next, if a is contracted with b, there 
are again two inequivalent possibilities for contracting c leading to (ab, ci, jk), (ab,ck,ij). Finally if a is contracted 
with k, there are three possible ways of contracting b giving (ak, bj, ci), (ak, be, ij), {ak, bi, cj). 

Of these, the contraction (ak, be, ij) is the same as (ab, ci,jk) since (ak, be, ij) = (ic, jk, ab) under (a, b, c) <-> (i,j, k) 
and, of course, (ic, jk,ab) — (ab,ci,jk). Similarly, (ak,bj,ci) — (ic,jb,ka) = (ib,jc,ka); the first equality comes 
from (a,b,c) <-> (i,j,k) symmetry while the second arises from b <-> c symmetry. Since (ib,jc,ka) = (ak,bi,cj) 
trivially, we need to retain only the first two out of the three possibilities in the last set. Thus dropping the 
two contractions (ab,ci,jk) and (ak,bi,cj) out of the 7 possibilities, we are left with 5 different contractions: 
(ai, bc,jk), (ai, bj, ck), (ab, ck, ij), (ak, bj, ci), (ak, be, ij). This will correspond to an action for the spin-2 field of the 
form: 

A=- I d i xd a h bc d i h jk [ciV ai V bc rf k + c 2 r 1 al r 1 b: >r 1 ck +c 3 v ab jf k ^+c 4 n ak v bl n cl +c 5 r l ak 7 1 bc j 1 13 ] 



d 4 x 



Cl d a h b b d a h] + c 2 d a h bc d a h bc + c^d a h ab d l W b + cM bc d c h ba + c 5 d a h b b dih ia (62) 



Each term in the action in Ea. (|62f> is of the kind d a hb c dihjkJ abcl ' : ' k (ri) where J abci 3 k is a cubic in rj and hence is 
constant [i.e, all components are or ±1]. This allows one to "swap" the derivatives di and d a by adding a total 
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divergence, using the identity: 

[d a h bc dih jk - dih bc d a h jk ]J abcijk = d a [h bc dih jk (J abcijk - r b ^ k )] 
Using this result, one can convert the C3 term to the C4 term and rewrite the action as 



(63) 



.4 



= - I d 4 x 



where 



^4div — 



d 4 xd a h bc dih jk [citfW* + c 2 v ai V bj V ck + (ca + c^^rj^ + c 5 r 1 ak r } bc r 1 ^] + A dn 



c x d a h b h d a h\ + c 2 d a h bc d a h bc + (ca + c 4 )d a h bc d c h ba + c 5 d a h b b d t h 



A div ee A h + A 



div 



(64) 



ab^ck^ij ^ak^bj '^ci] £3 

4 



d 4 xd a h bc dih jk [rfW - V ak rf^ 
j y d^xdaih^dihi-h^dih 



Jd 4 x [d a h ab d t hi-d a h bc d c h b 



(65) 



which, being a total divergence, does not contribute to the equations of motion if suitable boundary conditions are 
imposed. Notice that there are no further ambiguities related to "swapping" of derivatives in the action in Eq. I|62[) : 
this is clearly not possible in the c\,C2 or C5 terms, since the swapping leads to identical terms. Hence the only 
ambiguity is in the choice between the C3 term for C4 terms. 

Interestingly enough, the constants ci,C2,c 5 , (03 + C4) in Ah in Ea. (|64|l can be determined except for an overall 
scaling by the requirement that the field equations should be invariant under the gauge transformation: 



h ab (x) -> h ab (x) + d a £ b (x) + d b (, a (x). 



(66) 



This fixes the constants to be c\ = —c^ = 1; C3 + C4 = — C5 = 2 except for an overall scaling which is left as (1/4) for 
future convenience. The resulting expression for the quadratic part of the action can be written in different forms: 



A, 



d A xd a h bc dih jk [?7 a y V fe - r? a W fc + 2rf k r ] bj r ] cl - 2r} ak T] bc r} ij ] 



= iUx 



= A 



diK&Hi - d a h bc d a h bc + 2d a h bc d c h ba - 2d a h b d l h 



-3^8% - d a h bc d a h bc + 2d a h bc d c h ba 



h ab — h ab — -r] ab hl 



We shall use the shorter notation: 



A h = - t I d A xd a h bc dih jk M ab ^ k {r, mn ) 



(67) 



(68) 



where the tensor M abc ^ k (rj mn ) is symmetric in bc,jk and under the triple exchange (a,b,c) *-* (i,j,k) and is given 
by: 



symm 



(69) 



where the subscript "symm" indicates that the expression inside the square bracket should be suitably symmetrized 
in bc,jk and under the exchange (a, 6, c) «-> (i,j,k). In the expression for the action, since 

M a b djk ig mu i t ipii e d 

by d a h bc dihj k , we need not worry about symmetrization and use the expression given inside the square bracket in 
Eq. {§§1 as it is. 



B. Gauge conditions and true degrees of freedom 

The gauge invariance which was imposed to obtain the action in Eq. (|67|l implies that we are dealing with redundant 
degrees of freedom in h ab and — without additional restrictions — it does not carry pure spin-2, in the sense of 
irreducible representations of Lorentz group. To see this explicitly, consider a hi k of the form: 

h ik (x) = Q lk (x) + d l A k {x) + d k Ai{x) + (did k - \ mk dA a(x) + ^rj ik l3{x) (70) 
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where 

0iQ* = O; mkQ ik = 0; d t A l = (71) 

so that /life (10 components) is separated into a transverse traceless tensor (10 — 5 = 5 components), a transverse 
vector (4—1 = 3 components) and two scalars a and (3. The action in Eq. (|67[1 now becomes 



09 - d 2 a) (72) 



[This is most easily proved in the Fourier space using the third line of Eq. 1|67|) ]. This expression shows that: (a) 
the action is independent of the vector degree of freedom as one would have guessed from the gauge invariance; (b) 
it does depend on the scalar e(x) but e and Q ab are decoupled from each other; (c) the residual scalar appears with 
the wrong sign for the kinetic energy term. Therefore, isolating the physical degrees of freedom requires imposing the 
conditions: h a a — [to set e = 0] and d a h c ^ = [to ensure transverse-traceless condition on Q ab '.] If we impose these 
gauge conditions in the action itself, it becomes: 

Ah = ^J d 4 xd a h bc dih jk + 2r?Vy fc ] = \J^x [-d a h bc d a h bc + 2d a h bc d c h ba ) (73) 

We mentioned before that the original action in Eq. l|67|) there was an ambiguity with respect to the two terms C3 
and C4. This non uniqueness disappears in Ea. (|73|l . There is no possibility of "swapping" the derivatives in the first 
term [since both are d a }; in the second term, swapping the derivatives will give a vanishing term because d a h ba = 0. 
We shall see in Appendix B that this action generalizes to an interesting nonlinear theory. (The five degrees of 
freedom in Q ab can be further reduced to the two physical degrees of freedom of a graviton by using a residual gauge 
transformation of the form in Eq. jSBJ) with D£ a = 0; see e.g., page 946 of Q; we will not require this result in our 
discussion) . 



C. Gauge invariance and conservation of the source 



The symmetries of the theory are easier to see in the momentum space, which can be done by introducing the 
Fourier components f a b(p) of h a b(x) defined as usual by: 

K b {x) = J ^L/ a6 (p) e ** (74) 

The action becomes 

Ah = I / Wf h ^ [P^ MabCl3k (v mn )] - \ J ^-J b cf; k N bc > k (75) 

with 

N bcjk = j- p 2 (jjbc^jk _ ^yfc) + 2p k ^ pCrj bj _ pJrjbc^ 

= (pV V fe - pW c - p c pV*) - y (?7 b V fe + V CJ V bk ) + \ (p k p c V bj + P 3 p c V bk + P k p b V C3 + P 1 p b rf k ) (76) 

where we have exhibited the symmetrized expression in full glory for once. \N bc i k is symmetric in bc,jk and under 
the pair exchange (6,c) <-> (j, k).] The gauge transformation of the spin-2 field, given by Ea. (|66(l is equivalent to 
fab —> fab + Pa£,b + Pb£a in the Fourier space. Using this in Eq. (|75^l we get 

Ah \ / Wf Ubc + 2pb ^ f ^ + *P&) NbBjh ( ?? ) 

Straight forward computation now shows that p\,N bc ^ k = = pjN bc i k making Ah is invariant under the gauge 
transformations. This is, of course, built-in by the choice of the coefficients a in the original action. This condition 
translates, in coordinate space, to the identity M 3 dbd a dihjk = 0. 
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Appendix B: Comment on the uniqueness 

We briefly comment on two issues in this Appendix, postponing their detailed discussion to a future publication. 

First, it may be noted that the action for spin-2 graviton had an ambiguity in the form a four divergence term 
Adiv in Eq. I|65[l . While the action in Eq. (0J correctly leads to full Einstein theory under the substitution d A x — > 
d A x^f^g, rj a b — > g a i, = r] a b + Xh a b, the four divergence term in Eq.|j|2J does not map to a four divergence term under 
these substitutions: 

which cannot be expressed as a four divergence. So, we obtain the correct Einstein's theory only if we start with the 
correct set of terms in the linear spin-2 case. In general, consider two different actions Aj and A 1 ™ in the linear 
theory, which differ by a four divergence. The substitutions d 4 x — > d 4 Xy/—g, r\ a b — * 9ab — Vab + Xh a b can be used with 
either A l } n or A l }f and will lead to nonlinear theories with actions A™ 1 or A™\. But A 7 } 1 and A™\ cannot be related 
by a four-divergence, in general. 

Second, let us consider the spin-2 theory obtained by imposing the gauge conditions h® = 0, d a h^ — in the action 
itself. The linear theory action was given by Eq. I|73l) and, as was pointed out before, there is no ambiguity as regards 
four divergences in this action. Our "rule" now leads to the nonlinear action: 

A 00 = ±[d A xV=9 d a9bc d l9jk '/"//V" + 2 5 W1 (79) 

To get the field equations from this action, we only need to note that it is the same as the standard T 2 action in 
Eq Q23|l with the extra condition g a bdg ab — 0. Hence the field equations resulting from varying g a b in Eq. (|79|l will 
be the same as those obtained from varying g a b in the action in Eq l|23|l . keeping \J—g = constant. It is well known 
that, we then get the equations: 

Rab - jQabR = 87r(T ab - ^g a bT) (80) 

in which both sides are trace free. Bianchi identity can now be used to show that d a (R + 8irT) — 0, requiring 
(R + 8nT) — constant. Thus cosmological constant arises as an (undetermined) integration constant in such models 
fl4| . and could be interpreted as a Lagrange multiplier that maintains the condition y/—g = constant. Clearly, the 
gauge conditions translate to \J—g = constant in the full theory eliminating the scalar degree of freedom. 

It should be stressed that, while this theory is mathematically the same as Einstein's gravity with a cosmological 
constant, it is conceptually quite different from the usual approach to cosmological constant. What is probably more 
interesting, this theory takes gravity one step closer to gauge theories in the following sense: It has been known for a 
long time that the Christoffel symbols in gravity r^ b can be thought of as the elements of the matrix (r o )£ in exact 
analogy with the matrix representation of the gauge field (A^f^. The Riemann-Christoffel tensor can then be given 
the matrix representation 

Rab = d a T b - d b T a + T a T b - T a T b (81) 

with two matrix indices suppressed etc. Everyone who followed this route soon realized that — unfortunately — we 
need to contract on a matrix index with spacetime index to get the Einstein action etc. In the above approach, if we 
take \f—g = 1 then the quadratic action can be expressed in the form 

Auad - ^ / d4x 9 ab Tr[T a T b ] ; yf~g = 1 (82) 
We hope to discuss the gauge theory connection arising from this approach in a future publication. 
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